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ABSTRACT 

Low-luminosity active galactic nuclei (LLAGNs) represent the bulk of the AGN population in the 
present-day universe and they trace the low-level accreting supermassive black holes. In order to 
probe the accretion and jet physical properties in LLAGNs as a class, we model the broadband radio 
to X-rays spectral energy distributions (SEDs) of 21 LLAGNs in low-ionization nuclear emission-line 
regions (LINERs) with a coupled accretion-jet model. The accretion flow is modeled as an inner ADAF 
outside of which there is a truncated standard thin disk. We find that the radio emission is severely 
underpredicted by ADAF models and is explained by the relativistic jet. The origin of the X-ray 
radiation in most sources can be explained by three distinct scenarios: the X-rays can be dominated 
by emission from the ADAF, or the jet, or the X-rays can arise from a jet- ADAF combination in which 
both components contribute to the emission with similar importance. For 3 objects both the jet and 
ADAF fit equally well the X-ray spectrum and can be the dominant source of X-rays whereas for 1 1 
LLAGNs a jet-dominated model accounts better than the ADAF-dominated model for the data. The 
individual and average SED models that we computed can be useful for different studies of the nuclear 
emission of LLAGNs. From the model fits, we estimate important parameters of the central engine 
powering LLAGNs in LINERs, such as the mass accretion rate and the mass- loss rate in the jet and 
the jet power - relevant for studies of the kinetic feedback from jets. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galaxies: nuclei 
— galaxies: jets — galaxies: Seyfert 


1. INTRODUCTION 

One of the central paradigms in extragalactic astron- 
omy is that today’s galaxies host supermassive black 
holes (SMBHs) at their centers and have a symbiotic 
evolution with them (Magorrian et al. 1998; Ferrarese 
& Ford 2005). The central black holes seeded in proto- 
galaxies at high redshift grow in mass during cosmic 
history through a sequence of mergers of massive black 
holes and accretion episodes (e.g., Volonteri et al. 2005; 
Fanidakis et al. 2011). The massive black holes grow in 
such a way that they co-evolve with the galactic bulges 
(Kormendy k Richstone 1995; Ferrarese & Merritt 2000; 
Gebhardt et al. 2000; Ferrarese & Ford 2005) but do not 
seem to do so with the galaxy disks (Kormendy k Ben- 
der 2011) or dark matter haloes (Kormendy, Bender k 
Cornell 2011). 

The growih of SMBHs is dominated by the quasar 
phase (Marconi et al. 2004; McLure & Dunlop 2004; Si- 
jacki et al. 2007; Merloni k Heinz 2008) during w^hich 
they accrete from their large gas reservoir via geometri- 
cally thin accretion disks which are radiatively efficient 
(Shakura k Sunyaev 1973; Novikov & Thorne 1973) as 
suggested by different lines of evidence (Shields 1978; Ko- 
ratkar k Blaes 1999: Shang et al. 2005). The quasar 
population has a strong cosmological evolution with its 
density decreasing by a factor of ~ 100 from z ~ 2 - 
where the density reaches its peak - to the present time 
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(e.g., Osmer 2004; Silverman et al. 2005). 

In the present-day universe, SMBHs are underfed com- 
pared to the ones at high z and are “sleeping”. Most 
of SMBH activity at low z is dominated by the weak 
end of the AGN luminosity function in the form of low- 
luminosity AGNs (LLAGNs: Ho, Filippenko k Sargent 
1995, 1997b; Nagar, Falckefe Wilson 2005; Ho 2008), 
Though the LLAGN phase dominates the time evolu- 
tion of SMBHs, it contributes little to their mass growth 
(Hopkins k Hernquist 2006; Sijacki et al. 2007; Merloni 
k Heinz 2008; Xu k Gao 2010). 

The bulk of the LLAGN population (« 2/3; Ho 2008, 
2009) reside in low- ionization nuclear emission-line re- 
gions (LINERs; Heckman 1980; Ho, Filippenko k Sar- 
gent 1997a). LLAGNs are extremely sub-Eddington sys- 
tems which are many orders of magnitude less lumi- 
nous than quasars, with average bolometric luminosities 
(Lboi) ~ 10 40 ~~ 10 41 erg s”* 1 and an average Eddington 
ratio of Lboi/L^dd ~ 10” 5 (Ho 2009) where Lead is the 
Eddington luminosity. The observational properties of 
LLAGNs are quite different from those of more luminous 
AGNs. Regarding the broadband spectral energy distri- 
butions (SEDs), LLAGNs seem not to have the thermal 
continuum prominence in the ultraviolet (UV) - the “big 
blue bump” - which is one of the signatures of the pres- 
ence of an optically thick, geometrically thin accretion 
disk (Ho 1999; Nemmen et al. 2006; Wu, Yuan & Cao 
2007; Ho 2008; Eracleous, Hwang k Flohic 2010). Re- 
garding the emission-lines, LLAGNs typically have weak 
and narrow' Fe K a emission (Terashima et al. 2002) and 
a handful of LINERs display broad double-peaked Ha 
lines (e.g., Storchi-Bergmann et al. 2003); these proper- 
ties of the emission-line spectrum are consistent with the 
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absence of a thin accretion disk, or a thin accretion disk 
whose inner radius is truncated at > 100 GMf<? (Chen, 
Halpern k Filippenko 1989; Chen k Halpern 1989). Last 
but not least, with the typical fuel supply of hot diffuse 
gas (via Bondi accretion) and cold dense gas (via stellar 
mass loss) available in nearby galaxies, LLAGNs would 
be expected to produce much higher luminosities than 
observed on the assumption of standard thin disks with 
a 10% radiative efficiency (Ho 2009). Taken together, 
this set of observational properties favors the scenario 
in which the accretion flow in LLAGNs is advection- 
dominated or radiatively inefficient. 

Advection-dominated accretion flows (ADAFs 4 * & ; for re- 
views see Narayan, Mahadevan k Quataert 1998; Yuan 
2007; Narayan k McClintock 2008) are very hot, geo- 
metrically thick, optically thin flows which are typified 
by low radiative efficiencies (L <C O.lMc 2 ) and occur 
at low accretion rates (M < 0.01 Mead)* SMBHs are 
thought to spend most of their lives in the ADAF state 
(Hopkins & Hernquist 2006; Xu k Gao 2010), the best 
studied individual case being Sgr A* (e.g., Yuan 2007). 

In many LLAGNs, another component in the accretion 
flow besides the ADAF is required in order to account 
for different observations, including a “red bump’ 1 in the 
SEDs (Lasota et al. 1996; Quataert et al. 1999; Nem- 
men et ah 2006; Yu, Yuan k Ho 2011) and as mentioned 
before the double-peaked Balmer emission lines (e.g., 
Storchi-Bergmann et al. 2003; Eracleous, Lewis k Flohic 
2009): the emission from a thin accretion disk whose in- 
ner radius is truncated at the outer radius of the ADAF. 
In many LLAGNs the accretion flow may begin as a stan- 
dard thin disk but somehow at a certain transition ra- 
dius the accretion flow abruptly switches from a cold to 
a hot ADAF mode. The details of how this transition 
might happen are still not well understood (Manmoto 
et al. 2000; Yuan k Narayan 2004; Narayan k McClin- 
tock 2008), but it seems to be analogous to the transition 
between the different spectral states in black hole binary 
systems (Remillard k McClintock 2006; Done, Gierlinski 
k Kubota 2007). 

Maoz (2007) challenged the scenario of the central en- 
gines of LLAGNs consisting of ADAFs and truncated 
thin disks. Maoz argues that LLAGNs in LINERs have 
UV/X-rav luminosity ratios similar on average to those 
of brighter Seyfert Is and based on that observation he 
posits that thin disks extending all the way down to the 
radius of the innermost stable circular orbits (ISCO) per- 
sist even for LLAGNs, despite their smaller accretion 
rates. Yu, Yuan & Ho (2011) showed that the SEDs com- 
piled by Maoz (2007) are naturally fitted by ADAF mod- 
els with the addition of a truncated thin disk and also dis- 
cussed on theoretical grounds why the ADAF model has 
a superior explanation capability than the pure thin disk 
model. The ADAF model is the only model that can nat- 
urally account for the set of observational properties of 
LLAGNs within a self-contained theoretical framework. 
Hence, it is the physical scenario adopted in this work. 

ADAFs are relevant to the understanding of AGN feed- 
back since they are quite efficient at producing power- 

4 In this paper, we consider ADAFs and radiatively inefficient 

accretion flows (RIAFs) to be the same kind of accretion flow so- 
lution. For a clarification regarding the terminology, see Yu. Yuan 

& Ho 2011. 


ful outflows and jets, as suggested by theoretical stud- 
ies including analytical theory (Narayan k Yi 1994; 
Meier 2001; Nemmen et al. 2007) and numerical sim- 
ulations (McKinney & Gammie 2004; Hawley k Kro- 
lik 2006; Tchekhovskoy, Narayan k McKinney 2010)). 
This is in line with the different observational studies 
which demonstrate that LLAGNs are generally radio- 
loud (Ho 2002, 2008) and accompanied by powerful jets 
(e.g., Heinz, Merloni k Schwab 2007; Merloni k Heinz 
2008). In fact, the so-called “radio mode” of AGN feed- 
back invoked in semi-analytic and hydrodynamic simula- 
tions of galaxy formation (e.g., Bower et al. 2006; Croton 
et ah 2006; Sijacki et al. 2007; Fanidakis et al. 2011) 
would presumably correspond to the ADAF accretion 
state actively producing jets as explicitly incorporated 
in some works (Sijacki et al. 2007; Okamoto, Nemmen k 
Bower 2008). 

It is clear that an understanding of the physical na- 
ture of LLAGNs in LINERs will shed light on the nature 
of black hole accretion, outflows and consequently black 
hole growth and AGN feedback in present-day galaxies. 
One of the best ways of exploring the astrophysics of 
black hole accretion and outflows is by using multiwave- 
length observations of black hole systems and compar- 
ing the spectra predicted by specific models with the 
data. The goal of this work is therefore to probe the 
physics of accretion and ejection in the LLAGN popula- 
tion, by modeling their nuclear, broadband, radio to X- 
rays SEDs which provide constraints on physical models 
for the emission of the accretion flow and the jet. Fur- 
thermore, with a large enough sample of systems, we 
can derive from the fits to the data the parameters that 
characterize the central engines and build a census of the 
“astrophysical diet” of low-state AGNs. 

In this work, we will present how the currently fa- 
vored model for the central engines of weakly accret- 
ing black holes is able to account for the observed 
SEDs of LLAGNs. We will also present the typical ac- 
cretion/ejection - or feeding/feedback - properties of 
LLAGNs in LINERs as inferred from our detailed SED 
fits. With these parameters, we will then be able to 
draw a connection between the SMBHs in nearby active 
galaxies and the one in the Milky Way, Sgr A*. We will 
also present the inferred production site of the continuum 
emission and other properties inferred from the fits. 

The structure of this paper is as follows. In Section 2 
we describe the sample of 21 LLAGNs in LINERs that 
we used, the data and classifications. Section 3 describes 
the physical model that we adopted in order to interpret 
and fit the SEDs and derive the central engine parame- 
ters. In Section 4 we describe the model fits to the SEDs, 
explaining in more detail the SED models for 5 objects 
which we consider the “reference” or fiducial LLAGN ex- 
amples. Section 5 presents the model fits to the other 16 
SEDs in our sample. In Section 6 we describe the distri- 
butions of central engine parameters resulting from our 
detailed model fits. We present the average SED result- 
ing from our fits -- including the predicted emission in 
wavebands which can be observed with future facilities 
such as ALMA and the James Webb Space Telescope - 
and the inferred nature of the broadband emission in each 
wavelength in §7. We discuss in 8 which model is favored 
to explain our LLAGN sample, in particular which com- 
ponent in the flow dominates the X-ray emission. We 
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conclude by presenting a summary of our results in 9 
and comparing the SMBH “diet” of massive black holes 
in LINER AGNs and Sgr A*, 

2. SAMPLE 

Eracleous, Hwang & Flohic (2010) (hereafter EHF10) 
compiled a sample of 35 SEDs of LLAGNs found in LIN- 
ERs which include high spatial resolution optical and UV 
observations with the Hubble Space Telescope (HST), as 
well as X-ray observations with Chandra. Most of the 
SEDs studied by EHF10 have also high-resolution ra- 
dio (observed with the Very Large Array - VLA - or 
VLBA/VLBI) and infrared observations available. In or- 
der to select the best SEDs to be modelled from the sam- 
ple of EHF10, we applied the following selection criteria: 
we selected only those objects for which their X-ray flux 
is not an, upper limit. We also demand that there is a 
determination of the black hole mass. These selection cri- 
teria leave us with 24 LINERs. In three of these objects 
(NGC 404 , NGC 5055 and NGC 6500) the emission is 
likely to be dominated by a stellar population (EHF10). 
Hence we discard these three LLAGNs and we are left 
with 21 objects which are listed in Table 1. This table 
also lists their corresponding Hubble and LINER Types, 
black hole masses, bolometric luminosities and Edding- 
ton ratios. 

According to the quality of the available data, we can 
further classify our sample of LINERs in two groups. 
Group A contains the objects that have the data with 
the best quality available. For instance, the radio, op- 
tical and UV bands have each at least one data point 
which is not an upper limit to the luminosity. Obviously, 
the objects in this group provide better constraints to 
the accretion-jet models. Group B comprises the sources 
that provide not as good contraints to the models, be- 
cause they lack observational constraints in the radio to 
UV region of the spectrum. Some objects in this group 
have only upper limits to the flux in the UV, others have 
no data points or only upper limits in the radio band. All 
the objects in both groups have observations in the near- 
IR band which are treated as upper limits. We take this 
approach because the observations in the near-IR were 
taken with lower spatial resolutions, i.e. larger apertures 
that include considerable contamination from the emis- 
sion of the host galaxy. Therefore, near-IR observations 
can only be considered as upper limits to the nuclear 
emission. 

The masses of the SMBHs were estimated from the 
stellar velocity dispersions using the Mbh ~~ <? relation- 
ship (Ferrarese & Merritt 2000; Gebhardt et al. 2000; 
Tremaine et al. 2002), with the exception of NGC 266, 
whose black hole mass was estimated from the measured 
5 GHz radio and the 2-10 keV X-ray luminosities using 
the correlation of Merloni, Heinz & Di Matteo (2003) 
(the “fundamental plane of black hole activity”), and 
NGC 3031 (M81) and NGC 4486 (M87) whose black hole 
masses were estimated from the stellar and/or gas kine- 
matics (Bower et al. 2000; Devereux et al. 2003; Gebhardt 
et al. 2009; see Table 1). As EHF10 note, for the objects 
with multiple mass determinations using different meth- 
ods, the estimated masses are consistent with each other 
within a factor of 2. Due to the intrinsic scatter in the 
fundamental plane, the mass estimated for NGC 266 is 
subject to an uncertainty by a factor of & 10 (Merloni, 


Heinz & Di Matteo 2003). 

The optical-UV data of all objects were corrected for 
extinction as discussed in EHF10. In order to compute 
the bolometric luminosities from the SEDs, EHF10 used 
two methods. For the objects with the best sampled 
SEDs, they computed Lboi by integrating the SEDs di- 
rectly, ignoring upper limit data points. The objects 
with well-sampled SEDs are NGC 3031, NGC 3998, NGC 
4374, NGC 4486, NGC 4579 and NGC 4594. For these 
objects they assumed that pairs of neighboring points in 
the SEDs defined a power law, integrated each segment 
individually and summed the segments to obtain Lboh 
From this set of best sampled SEDs they estimated the 
average “bolometric correction” from the 2-10 keV lu- 
minosity to the bolometric luminosity (Lboi = 50 Lx), 
which they used to obtain Lboi for the remaining ob- 
jects, for which the SEDs are not as well sampled. These 
bolometric luminosities are listed in Table 1. 

3. MODELS FOR THE ACCRETION FLOW AND 
JET 

We fit the observed broadband SEDs of the LLAGNs in 
our sample using a model which consists of three compo- 
nents: an inner ADAF, an outer truncated thin accretion 
disk and a jet. The components of the model are illus- 
trated in Figure 1. We describe here the main features 
of this model. 



Figure 1 . Cartoon illustrating the model for the central engines 
of LLAGNs. It consists of three components: an inner ADAF, an 
outer truncated thin disk and a relativistic jet. 


3.1. ADAF component 

The inner part of the accretion flow is in the form of an 
ADAF which is a hot, geometrically thick, optically thin 
two- temperature accretion flow, which has low radiative 
efficiency (e.g., Narayan, Mahadevan & Quataert 1998; 
Kato, Fukue & Mineshige 1998). ADAFs are character- 
ized by the presence of outflows or winds, which prevent 
a considerable fraction of the gas that is available at large 
radii from being accreted onto the black hole. This has 
been suggested by numerical simulations (Stone et al. 
1999: Hawley, Balbus & Stone 2001; Igumenshchev et 
al 2003; De Villiers, Hawley & Krolik 2003; McKinney 
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Table 1 

Sample of galaxies and their basic properties® 


Galaxy 

Hubble 

Distance 15 

log 

LINER 

Lx 

L b oi 

^bol/^Edd SED 


Type 

(Mpc) 

(Mbh/Mq) 

Type 

(ergs r ) c 

(erg s “ 1 ) d 

type® 


NGC 266 

SB(rs)ab 

62.4 

(i) 

7.6 f 

LI 

7.4 x 

10 40 

2.2 x 10 42 

4 x 10“ 4 

B 

NGC 1097 

SB(rl)b 

14.5 

a) 

8.1 

LI 

4.3 x 

10 40 

8.5 x 10 41 

5 x 10"' 5 

A 

NGC 1553 

SA(rl)0 

17.2 

(2) 

7.9 

L2/T2 

1.3 x 

10 40 

3.8 x 10 41 

4 x 10" 5 

B 

NGC 2681 

SBA(rs)0/a 

16.0 

(2) 

7.1 

LI 

6.1 x 

10 38 

1.8 x 10 40 

1 x 10- 5 

B 

NGC 3031 (M81) 

SA(s)ab 

3.6 

(3) 

7.8 

SI. 5/Ll 

1.9 x 

10 40 

2.1 x 10 41 

3 x 10" 5 

A 

NGC 3169 

SA(s)a 

19.7 

(1) 

7.8 

L2 

1.1 x 

10 41 

3.3 x 10 42 

4 x 10~ 4 

B 

NGC 3226 

E2 

21.9 

(2) 

8.1 

LI 

5.0 x 

10 40 

1.5 x 10 42 

1 x 10”“ 4 

B 

NGC 3379 (M105) 

El 

9.8 

(2) 

8.2 

L2/T2 

1.7 x 

10 37 

5.1 x 10 38 

3 x 10~ 8 

B 

NGC 3998 

SA(r)0 

13.1 

(2) 

8.9 

LI 

2.6 x 

10 41 

1.4 x 10 43 

1 x 10~ 4 

A 

NGC 4143 

SAB(s)0 

14.8 

(2) 

8.3 

LI 

1.1 x 

10 40 

3.2 x 10 41 

1 x 10~ 5 

A 

NGC 4261 

E2-3 

31.6 

(2) 

8.7 

L2 

1.0 x 

10 41 

6.8 x 10 41 

1 x 10~~ 5 

B 

NGC 4278 

El-2 

14.9 

(2) 

8.6 

LI 

9.1 x 

10 39 

2.7 x 10 41 

5 x 10~ 6 

A 

NGC 4374 (M84, 3C 272.1) 

El 

17.1 

(2) 

8.9 

L2 

3.5 x 

10 39 

5.0 x 10 41 

5 x 10~ 6 

A 

NGC 4457 

SAB(s)0/a 

10.7 

(4) 

6.9 

L2 

1.0 x 

10 39 

3.0 x 10 40 

3 x 10" 5 

B 

NGC 4486 (M87, 3C 274) 

E0-1 

14.9 

(2) 

9.8 

L2 

1.6 x 

10 4 ° 

9.8 x 10 41 

7 x 10~ 6 

A 

NGC 4494 

El-2 

15.8 

(2) 

7.6 

L2 

9.2 x 

10 38 

2.8 x 10 4 ° 

6 x 10”~ 6 

B 

NGC 4548 (M91) 

SBb(rs) 

15.0 

(3) 

7.6 

L2 

5.4 x 

10 39 

1.6 x 10 41 

3 x 10~ 5 

B 

NGC 4552 (M89) 

E 

14.3 

(2) 

8.2 

T2 

2.6 x 

10 39 

7.8 x 10 4 ° 

4 x 10“ 6 

A 

NGC 4579 (M58) 

SAB(rs)b 

21.0 

(4) 

7.8 

L2 

1.8 x 

10 41 

1.0 x 10 42 

1 x 10~ 4 

A 

NGC 4594 (M104) 

SA(s)a 

9.1 

(2) 

8.5 

L2 

1.6 x 

10 4 ° 

4.8 x 10 41 

1 x 10~ 5 

A 

NGC 4736 (M94) 

(R)SA(r)ab 

4.8 

(2) 

7.1 

L2 

5.9 x 

10 38 

1.8 x 10 40 

1 x 10" 5 

A 


a The information in this table was obtained from EHF10, see text. 

k The number in parenthesis gives the source and method of the distance measurement, as follows: (1) From the catalog of Tully (1988), 
determined from a model for peculiar velocities and assuming Ho — 75 km s'* 1 Mpc" 1 ; (2) From Tonry et al. (2001), who used the surface 
brightness fluetation method. Following Jensen et al. (2003), the distance modulus reported by Tonry et al. (2001) was corrected by 
subtracting 0.16 mag; (3) From Freedman et al. (1994, 2001), who used Cepheid variables; (4) From Gavazzi et al. (2009) who used the 
TuI ly- Fisher method. 

c Lx is the X-ray luminosity in the 2 — 10 keV range. 

d The bolometric luminosities of NGC 1097, NGC 3031, NGC 3998, NGC 4374, NGC 4486, NGC 4579 and NGC 4594 were estimated by 
integrating the SEDs. For all other galaxies Lboi was determined by scaling Lx as described in §2. 
e Classification of the SED depending on the quality of the data sampling (see Section 2). 

^ The mass estimate for NGC 266 is subject to a large uncertainty (fe 1 dex) since it is based on using the fundamental plane of Merloni, 
Heinz & Di Matteo (2003) (see text). 


& Gammie 2004) and analytical work (Narayan & Yi 
1994; Biandford & Begelman 1999; Narayan et al. 2000; 
Quataert & Gruzinov 2000a). In addition, observational 
studies suggest that winds coming from the inner parts 
of accretion flows are ubiquitous in AGNs (e.g., Tombesi 
et al. 2010a,b) and it is natural to take them into account 
when studying the LLAGN population. 

In order to take this mass- loss into account, we follow 
Biandford & Begelman (1999) and introduce the param- 
eter s by 

M - Mont ( TT~) , (1) 

to describe the radial variation of the accretion rate M out 
measured at the outer radius of the ADAF, R 0 ut . Follow- 
ing the results of numerical simulations of the dynamics 
of ADAFs, in our models we allow s to vary over the 
range 0 < s < 1, taking into account the current uncer- 
tainty in the determination of the amount of gas that is 
lost via winds in the accretion flow. 

The other parameters that describe the ADAF solu- 
tion are the black hole mass M ; the viscosity parameter 
a: the modified plasma 5 parameter, defined as the ratio 
between the gas and total pressures, 8 = P g /Ptou S, the 
fraction of energy dissipated via turbulence that directly 
heats electrons; and the adiabatic index 7. In our cal- 
culations, we adopt a = 0 . 3 , ,5 — 0.9 and 7 — 1.5. Tra- 


ditional ADAF models adopted <5 to be small (5 < 0.01; 
e.g., Narayan & Yi 1995). On the other hand, it has been 
argued that the value of S can be potentially increased 
due to different physical processes - such as magnetic re- 
connection - that affect the heating of protons and elec- 
trons in hot plasmas (e.g., Quataert & Gruzinov 1999; 
Sharma et al. 2007). Given the theoretical uncertainty 
related to the value of d, we allow it to vary over the 
range 0.01 < 5 < 0.5. 

The case in which <5 — 0.5 and s — 0.3 deserves spe- 
cial attention, since these values were obtained from fit- 
ting current ADAF models to the SED of Sgr A* (Yuan, 
Quataert & Narayan 2003; Yuan, Shen & Zhi-Qiang 
2006). We note that while Yuan, Yu & Ho (2009) fix 
the values of S and s in their fits, we leave these parame- 
ters free taking into account the theoretical uncertainties. 
There are degeneracies between s and S when fitting the 
data with ADAF models (Quataert et al. 1999), such that 
different combinations of the values of these parameters 
may result in similar SEDs. 

The cooling mechanisms incorporated in the calcula- 
tions are synchrotron emission, bremsstrahlung and in- 
verse Comptonization of the seed photons produced by 
the first two radiative processes. Given the values of 
the parameters of the ADAF, in order to compute its 
spectrum we first numerically solve for the global struc- 
ture and dynamics of the flow, as outlined in Yuan et ai. 
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(2000); Yuan, Quataert & Narayan (2003). Obtaining 
the global solution of the differential equations for the 
structure of the accretion flow is a two-point boundary 
value problem. This problem is solved numerically using 
the shooting method, by varying the eigenvalue j (the 
specific angular momentum of the flow at the horizon) 
until the sonic point condition at the sonic radius R s is 
satisfied, in addition to the outer boundary conditions. 

There are three outer boundary conditions that the 
ADAF solution must satisfy, specified in terms of the 
three variables of the problem: the ion temperature XJ , 
the electron temperature T e and the radial velocity v (or 
equivalently the angular velocity O). Following Yuan, 
Ma & Narayan (2008), when the outer boundary of the 
ADAF is at the radius R 0 ut ~ 1 0 4 Rs (where Rs is 
the Schwarzschild radius) we adopt the outer bound- 
ary conditions T ou u = 0.2T V i r , T ou t,e = 0.19T v i r and 
A ou t = 0.2, where the virial temperature is given by 
T v j r = 3.6 x 10 12 (it*s/i?) K, A ~ v/c s is the Mach 
number and c s is the adiabatic sound speed. When 
the outer boundary is at R ou t ^ 10 2 Rs we adopt the 
boundary conditions T out ,i = 0.6T v i r , T out , e — 0.08T V i r 
and A out = 0.5. After the global solution is calculated, 
the spectrum of the accretion flow is obtained (see e.g., 
Yuan, Quataert & Narayan 2003 for more details). We 
verified that if these boundary conditions are varied by 
a factor of a few, the resulting spectrum does not change 
much. 

3.2. Thin disk component 

Outside the ADAF there is an outer thin accretion 
disk with an inner radius truncated at R tr = R 0 u t and 
extending up to 10 5 i?s. Therefore, in practice the outer 
radius of the ADAF corresponds to the transition radius 
to the thin disk. The other parameters that describe 
the thin disk solution are the inclination angle i, the 
black hole mass and the accretion rate M out (the same as 
the accretion rate at the outer boundary of the ADAF). 
In the cases where R 0 ut ~ 10 A Rs w r e simply ignore the 
contribution of the thin disk spectrum. 

The thin disk emits locally as a blackbody and we take 
into account the reprocessing of the X-ray radiation from 
the ADAF. This reprocessing effect has only a little im- 
pact on the spectrum of the thin disk though, with the 
resulting SED being almost identical to that of a stan- 
dard thin disk (e.g., Frank, King & Raine 2002). 

3.3. Jet component 

The SEDs of LLAGNs are generally radio- loud (Ho 
1999; Ho & Peng 2001; Ho 2002; Terashima & Wilson 
2003; EHF10); but see Maoz 2007; Sikora, Stawarz & 
Lasota 2007). The radio prominence of these SEDs is 
usually explained by invoking the synchrotron emission 
of relativistic jets, since the the accretion Bow does not 
produce enough radio emission to account for the ob- 
served radio luminosity (e.g., Quataert et al. 1999; Ui- 
vestad & Ho 2001; Wu & Cao 2005; Nemmen et al. 2006; 
Wu, Yuan & Cao 2007; Yuan, Yu & Ho 2009). Some au- 
thors even argue that the entire SED of LLAGNs may be 
explained by the jet component (e.g., Falcke, Kording & 
Markoff 2004; Markoff et al. 2008). We therefore include 
in our modelling the contribution from the emission from 
a relativistic jet. 


We adopt a jet model based on the internal shock sce- 
nario which is used to interpret gamma-ray burst after- 
glows (e.g., Piran 1999; Spada et al. 2001; see Yuan, Cui 
& Narayan 2005 for more details). This model has been 
adopted in previous works to understand the broadband 
SEDs of X-ray binaries and AGNs (Spada et al. 2001; 
Yuan, Cui & Narayan 2005; Nemmen et al. 2006; Wu, 
Yuan & Cao 2007). According to this jet model, some 
fraction of the material in the innermost regions of the 
accretion flow is transferred to the jet producing an out- 
flow rate Mj e t and a standing shock wave in the region 
of the jet closest to the black hole is formed. This shock 
wave is created from the bending of the supersonic ac- 
cretion flow near the black hole in the direction of the 
jet. We calculate the shock jump (Rankine-Hugoniot) 
conditions to find the electron and ion temperatures of 
the plasma ( T e and Ti). We find that the jet spectrum 
is not very sensitive to changes in T e and T*, since the 
emission is completely dominated by non-thermal elec- 
trons (see below). We therefore adopt Ti = 6.3 x 10 11 K 
and T e — 10 9 K in our calculations of the jet emission. 

The jet is modelled as having a conical geometry with 
half-opening angle 4> and a bulk Lorentz factor T^ which 
are independent of the distance from the black hole. The 
jet is along the axis of the ADAF and makes an angle i 
with the line of sight. The internal shocks in the jet are 
presumably created by the collisions of shells of plasma 
with different velocities. These shocks accelerate a small 
fraction £ e of the electrons into a power-law energy dis- 
tribution with index p . The energy density of accelerated 
electrons and the amplified magnetic field are described 
by two free parameters, e e and e#. Following Nemmen 
et al. (2006); Wu, Yuan & Cao (2007), we adopt in our 
calculations of the jet emission the values <p = 0.1 ra- 
dians, = 10% and Fj = 2.3, which corresponds to 
v/c « 0.9 (except in the case of M87, which has an in- 
dependent estimate of Tj available). Unless otherwise 
noted, we adopt i = 30°. Therefore, there are four free 
parameters in the jet model: Mj e t? P, t e and cb* In our 
calculations we consider the synchrotron emission of the 
jet, with the optically thick part of the synchrotron spec- 
trum contributing mainly in the radio and the optically 
thin part contributing mainly in the X-rays. 

3.4. Free parameters and fit procedure 

When fitting the observed SEDs with our coupled 
accretion-jet model, we have eight free parameters. Four 
of these parameters describe the emission of the accretion 
flow: the accretion rate M out , the transition radius be- 
tween the inner ADAF and the outer truncated thin disk 
i? tr , the fraction of viscously dissipated energy that di- 
rectly heats the electrons S and the strength of the wind 
from the ADAF s. The other four parameters character- 
ize the jet emission: the mass-loss rate in the jet Mj e t, 
the electron energy spectral index p, and the electron and 
magnetic energy parameters e € and e^* M out and Mj et 
are obviously expected to be correlated, but given our 
ignorance about the mechanism of jet formation we vary 
these parameters independently in our fits. We require 
for consistency that Mj et /M ollt < 1. Since the ADAF 
loses mass via winds (equation 1), it is also of interest to 
check the value of Afi et /M(SRs), in which the accretion 
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rate is calculated at the radius of the innermost stable 
circular orbit for a Schwarzschild black hole. 

Throughout this paper we will use the dimension- 
less mass accretion rates m = M/Me dd? noting that 
the Eddington accretion rate is defined as Me <m = 
22M / (10 9 Mq) Mqjt^ 1 assuming a 10% radiative effi- 
ciency. We will also express the black hole mass in terms 
of the mass of the sun, m = M/M@, and the radius in 
terms of the Schwarzschild radius, r = R/R $ . 

Our fitting procedure can be summarized as follows. 
We vary the free jet parameters in order to fit the radio 
observations, since the radio is emitted predominantly by 
the jet. Once we fit the radio band with the jet model, 
we try to fit the truncated thin disk model to the IR- 
optical data, in order to constrain the transition radius 
and M out (keeping in mind that the theoretical spectrum 
should not exceed the observed luminosities). We then 
use the X-ray data to constrain the ADAF component 
in order to refine the estimates of R tr and M out , and 
to estimate s and S . Depending on the combination of 
parameters that we use, the jet can have an important 
contribution to the emission not only in the radio, but 
also in the X-ray band. We sum the emission of the 
inner ADAF, outer thin disk and jet, and compare this 
sum with the observed SEDs. 

We explore the parameter space of the models consid- 
ering the many possible fits to the data. We search the 
literature for independent constraints from other works 
on the values of the model parameters, such as the mass 
accretion rate, inclination angle and the transition ra- 
dius, as well as the jet power which is obtained from the 
jet model. 

4. MODEL FITS TO WELL SAMPLED SEDS 

We describe in this section the results obtained from 
fitting our coupled accretion-jet model to the SEDs that 
were selected using our selection criteria in Section 2. We 
describe in more detail the results for NGC 4594 (§4.1), 
NGC 4374 (§4.2), NGC 4486 (§4.3), NGC 3031 (§4.4) 
and NGC 3998 (§4.5). These five objects will serve as 
the “reference” examples that will illustrate the general 
results for the whole sample. These objects were cho- 
sen for a more detailed discussion because of their good 
multiwavelength sampling. We show the model fits to 
the other 16 SEDs and the resulting model parameters 
in Section 5. 

In the SED plots that follow below, the error bars rep- 
resent the uncertainty in the extinction corrections ap- 
plied to the the optical-UV data by EHF10. The lower 
bars correspond to the measurements without any ex- 
tinction correction, while the upper bars represent the 
same observations after a maximal extinction correction 
(see EHF10), The arrows represent upper limits to the 
nuclear luminosity. These upper limits are caused by 
either non-detections, or because the corresponding ob- 
servations were taken at lower spatial resolutions. In 
the latter case, the upper limits include a potentially 
significant contamination by the emission from the host 
galaxy. This happens particularly for the observations in 
the near-IR band. 

Table 2 lists the model parameters that result from 
the SED fits for the AD and JD type of models. We 
list whenever available the Bondi accretion rate and the 


corresponding reference from which meondi was taken. 

corresponds to the jet power estimated either from 
observations or using the correlation between jet power 
and radio luminosity of Merloni & Heinz (2007). Fj“ od 
represents the jet power resulting from our jet model, 
calculated as Pj e t = Mj et F 2 c 2 . 

4.1. The case of NGC 4594 

Pellegrini (2005) studied Chandra X-ray observations 
of NGC 4594 and derived the Bondi accretion rate 
TOBondi = 2±?;| X 10- 3 . 

The SED of NGC 4594 is plotted in Figure 2, together 
with two different spectral fits. The left panel of Figure 
2 shows an ADAF-jet. model in which the ADAF com- 
ponent (dashed line) completely dominates the emission 
from the IR to the X-ray bands and is consistent with the 
available optical-UV and X-ray data. In this model, the 
jet ( dot-dashed line) dominates the radio emission and 
contributes only weakly to the X-ray flux. The solid line 
corresponds to the sum of the ADAF and jet emission. 

The parameters of the ADAF model are m ou t = 
6.3 x io~ 3 (consistent with the Bondi accretion rate es- 
timated by Pellegrini 2005), r out = 10 4 , 6 = 0.01 and 
s = 0.3. The jet parameters are rhj e t = 9 x 10~ 7 , p = 2.3, 
e e = 5 x 10” 3 and e# = 0.03. The jet kinetic power is 
jFj e t/Lbol ~ 4. In other words, the kinetic power carried 
by particles in the jet exceeds the overall radiative power. 
The ratio of the mass-loss rate in the jet to the accretion 
rate estimated at 3 Rs is m^ et /rh(3Rs) ~ 1.6 x 10“ 3 . In 
other words, only « 0.2% of the mass that is ultimately 
accreted by the black hole is channeled into the jet. 

The dotted line in the left panel of Figure 2 shows an 
ADAF + truncated thin disk model for which r tr = 300, 
rhout = 9.1 x 10 -4 , S = 0.01 and s = 0.1. The thin 
disk emission corresponds to the “bump” at 10/im. This 
model demonstrates that accretion flows with a much 
smaller transition radius are also consistent with the 
data. Since we have only upper limit points in the mid 
to near-IR for NGC 4594, the available data are not suf- 
ficient to constrain the presence of a truncated thin disk 
in this source. 

The right panel of Fig. 2 shows a model for the SED 
of NGC 4594 in which the jet dominates the radio and 
X-ray emission. In fact, the right panel shows that the 
relativistic jet can account well for the entire observed 
SED . The jet parameters are mj et — 4.5 x 10“ 7 , p = 2.01, 
e e — 0.8 and e# = 3xl0 -2 , with Pj et = 8.4 xlO 42 erg s”” 1 , 
Pjet/^boi « 4 and m jet /m(3Ps) » 0.009. 

We display in the right panel of Figure 2 one illustra- 
tive ADAF model which is consistent with the possibility 
that the jet emission dominates in this source, i.e. the 
ADAF contribution in all bands (except at 1 mm) is very 
small compared to the jet. The ADAF parameters are 
m ou t “ hz Bondi 5 Tout ~ 10 4 , d — 0.01 and s = 0.3. We can 
see that the ADAF contribution in this case is only sig- 
nificant in the wavelength ranges ~ 1 mm — 100 pm. The 
importance of the ADAF emission in the other bands is 
small compared to the jet. 

By exploring the parameter space of the accretion /jet 
models within the range of plausible values allowed by 
theory, we demonstrate that there are basically two pos- 
sible types of models which can accommodate the ob- 
served SED of NGC 4594 and the other LINERs in 
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Figure 2. Models for the SED of NGC 4594 (Sombrero). The dashed and dot-dashed lines show respectively the emission from the 
ADAF and jet, while the solid line shows the sum of the radiation from these components. Left: model in which the ADAF dominates 
the observed X-ray emission. The dotted line shows for comparison an ADAF + thin disk model in which the transition radius is 300 Rs- 
Right: model in which the jet dominates the X-ray output. See the text for the parameters. 


this section. In the first type, the emission from the 
ADAF dominates the observed X-rays; in the second type 
of model, the jet emission dominates the X-rays. We 
will hereafter use the abbreviations AD (as in ADAF- 
dominated) and JD (as in jet-dominated) when referring 
to the former and latter types of models, respectively. A 
third type of model is also possible in which the jet and 
the ADAF contribute with similar intensities to the high 
energy emission. These results apply not only to the five 
LINERs discussed in more detail in this section but also 
to the other sources in our sample with a few exceptions. 

4.2. NGC 4374 

A prominent jet resolved with the VLA is observed 
to create cavities in the X-ray emitting gas (Allen et al. 
2006; Finoguenov et al. 2008). The Bondi accretion rate 
is mBondi = 4 x 10~~ 4 (Pellegrini 2005; Allen et al. 2006). 

The left panel of Figure 3 shows an ADAF- jet model 
( dashed line) in which the ADAF completely dominates 
the emission from the near-IR to the X-ray bands and 
explains well the available optical-UV and X-ray data, 
as was the case for NGC 4594. The synchrotron peak 
of the ADAF also dominates the emission in the range 
^ 1 mm — 100 jam. The jet ( dot-dashed line) dominates 
the radio emission, but contributes weakly to the X-ray 
flux. The thin disk dominates the IR emission and its 
contribution peaks at A ^ 10 ^m, creating a “red bump” 
in the SED as opposed to the “big blue bump” observed 
in quasars (e.g\, Koratkar & Blaes 1999; Nemmen & 
Brotherton 2010). The parameters of the accretion flow 
are m out = 3.95 x 10“ 4 , r tr = 150, 5 — 0.01 and s = 0.1. 
As was the case for NGC 4594, the mid to near-IR data 
are insufficient to make the case of a truncated thin disk 
compelling and similar models with much larger values 
of r tr are not ruled out by the available IR data. 

The jet parameters are m j et — 4 x 10~ 7 , p ~ 2.4, e e — 
0.01 and es = 0.1. The resulting jet power is Pj e t = 2.1 x 
10 42 erg s” 1 « 4Lhoi which is in excellent agreement with 


the jet power estimated by Merloni & Heinz (2007) based 
on the calorimetry of the X-ray cavities observed with 
Chandra (3.9 x 10 42 erg s'" 1 ). The ratio rh- ]e t/m(3Rs) is 
1.5 x 10" 3 . 

The right panel of Fig 3 shows a model in which the 
jet dominates the radio and the 1-100 keV X-ray band. 
The emission in the region between 1000 A and 1 keV is 
dominated by the Compton peak of the ADAF spectrum. 
As was the case for the previous model, the IR emission 
is dominated by a truncated thin disk with a small tran- 
sition radius (r tr = 30) which is even smaller than the 
AD fit. At v > 100 keV the bremmstrahlung emission of 
the ADAF is the dominant radiation process. The pa- 
rameters of ADAF are m ou t = 1.5 x 10“ 4 (in agreement 
with the Bondi rate) and S = s = 0.3. The jet parame- 
ters are m- ]et = 1.6 x 10“ 6 , p = 2.2, e e = 9 x 10“ 3 and 
cb = 8x 10~ 3 , with Pj e t = 8.4 x 10 42 erg s” 1 « 17Pboi* 
The ratio mj et /m(3Ps) is mj et /m(3Ps) « 0.02. There 
are no data available to constrain the predicted ADAF 
bumps at soft and hard X-rays in this model. Notice 
that, similarly to the JD fit to the SED of NGC 4594, the 
shape of X-ray spectrum from the ADAF wnth 5 = 0.3 
(dashed line) does not agree with the data. 

4.3. NGC 4486 

Di Matteo et al. (2003) previously calculated the Bondi 
accretion rate using Chandra X-ray data and found it 

to be MBondi ~ o.l M 0 yr -1 (m Bo ndi « 7 x 1CT 4 ). 
They also fitted the multiwavelength SED of M87 with 
an ADAF model for different values of S (but not in- 
cluding mass-loss in the accretion flow, i.e. s = 0) 
and found that the ADAF emission roughly explains the 
SED. Biretta, Sparks & Macchetto (1999) analysed HST 
observations of the jet in M87 and estimated that Tj > 6 
and 10° < i < 19°. Based on the results of Biretta, 
Sparks & Macchetto (1999), w^e adopt in our jet mod- 
elling the parameters Tj = 6 and i = 10°. The kinetic 
power carried by the jet is estimated to be in the range 
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Figure 3. Models for the SED of NGC 4374/M84 showing the emission of the ADAF ( dashed ), jet ( dot-dashed ), truncated thin disk 
( dotted ) and the total emission (solid). Left: model in which the ADAF dominates the observed X-ray emission (“AD model”). Right: 
model in which the jet dominates the X-ray output (“JD model”). 


10 43 — 10 44 erg s” 1 (e.g,, Bicknell & Begelman 1999; 
Owen, Eilek & Kassim 2000; Allen et al. 2006; Merloni 
& Heinz 2007), therefore the jet power resulting from 
fitting our jet model must be within this range. 

The left panel of Figure 4 shows an ADAF-jet model 
in which the ADAF dominates the optical-UV (OUV) 
and X-ray emission (i.e. AD-type model), while the jet 
dominates the radio band. The parameters of the ac- 
cretion flow are m out = 5.5 x 10“ 4 (in agreement with 
the Bondi rate), r out = 10 4 , 5 = 0.01 and s = 0.1. 
The jet parameters are rhj et = 5 x 10" 8 , p = 2.6, 
e e = 10~~ 3 and = 8 x 10~ 3 . The jet power is given 
by Pjet = 6 x 10 42 erg s” 1 , with Pj e t/Lboi ~ 6 and 
mj et /m(3i?s) «9x 10“ 5 . 

The right panel in Fig 4 shows a jet model which ex- 
plains quite well the radio and X-ray observations (JD- 
type model), although it underpredicts the OUV data 
by a factor of a few. The jet parameters are mj e t — 
3 x 10“~ 8 , p = 2.3, e e = 10” 3 and = 8 x 10“ 3 , with 
jfjet = 9 x 10 42 erg s _1 and P je t/iboi ~ 9. 

The SED of M87 was previously fitted using ADAF 
and/or jet models (Di Matteo et al. 2003; Yuan, Yu & Ho 
2009; Li et al. 2009). Di Matteo et al. (2003) modelled the 
SED of M87 with an ADAF model using different values 
of S but not including mass-loss (i.e., s — 0). The model 
adopted by Li et al. (2009) is quite similar to that of 
Di Matteo et al. (2003) although the former incorporate 
general relativistic corrections. Di Matteo et al. (2003); 
Li et al. (2009) obtained that the ADAF emission with no 
mass-loss approximately reproduces the SED and results 
in accretion rates consistent with the Bondi rate. Our 
AD model is similar to that of Di Matteo et al. (2003) 
but it also incorporates the effect of winds, as suggested 
by numerical simulations. 

Yuan, Yu & Ho (2009) tried to model the SED using 
an ADAF model with 5 = 0.5 but failed to fit to data 
with an AD model. The reason is that for high values of 
5 (S > 0.1) and small accretion rates (m out 0.01) the 
ADAF X-ray spectrum is harder than the data. Yuan, 


Yu & Ho (2009) instead successfully fit the data with a 
JD model using p = 0.5. 

4.4. NGC 3031 

This source is among the brightest LLAGNs known 
since it is the nearest AGN besides Centaurus A and 
has been the subject of a broadband multi wavelength 
monitoring campaign (Markoff et al. 2008; Miller et al. 
2010). Due to the many available observational con- 
straints, accretion-jet models can be well constrained for 
this LLAGN. For illustrative purposes, here we will re- 
port our modelling results using the data compiled by 
EHF10. 

Devereux & Shearer (2007) modeled the profile of the 
broad double-peaked Ha line with a relativistic thin disk 
model. They were able to explain the profile with an 
inclination angle of 50°. Such a high inclination angle is 
supported by radio observations of the jet (Bietenholtz, 
Bartel & Rupen 2000). The Ha line profile of M81 is 
consistent with an inner radius of « 280 — 360 Rs for the 
line-emitting thin disk (Devereux & Shearer 2007). In 
our models for the SED. we therefore adopt an inclination 
angle i = 50° and we impose the constraint on the models 
that r tr < 360. 

The left panel of Figure 5 show-s the AD fit model for 
the SED of M81. In this model, the synchrotron radi- 
ation from the jet is the dominant mechanism of radio 
emission. In the wavelength range between ~ 10 pm and 
« 1 pm a u red bump" can be seen which corresponds to 
the emission of the truncated thin disk. The thin disk 
dominates the emission in this band. In the interval be- 
tween 1 mm and ^100 pm the ADAF is the dominant 
component. The ADAF Is also the flow component that 
dominates the emission for wavelengths less than 1 p un- 
til X-rays. The sum of the radiation of all components is 
plotted as the solid line. The accretion flow parameters 
are m ou t = 3 x ICT 3 , r tr = 360 (consistent with the Ha 
modeling results), 5 = 0.01 and s = 0.16. The jet param- 
eters are m- }et == 2x 10~“ 6 , p — 2.2, e e = 0.1 and e# = 0.01. 
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The jet power is given by Pj et = 8x 10 41 erg s 3 with 
Pjet/Pbol « 4 and m^ et /m(SRs) « 10" 3 - 

In this AD model, a low value of S = 0.01 is required in 
order to reproduce the X-ray spectral shape character- 
ized by T = 1.88. As was the case for M87, if we increase 
the value to <5 = 0.3 keeping the other parameters fixed 
except s = 0.4 the predicted value of the photon index is 
reduced and the ADAF spectrum becomes harder than 
the data. 

The right panel of Fig, 5 shows a JD model in 
which the jet component dominates the radio, optical- 
UV and X-ray bands of the SED. In the wavelength 
range between 1 mm and « 100 pm the ADAF is the 
dominant component, although there are no observa- 
tional constraints to the model in this region. Between 
~ 10 pm and 6000 A the truncated thin disk dominates 
the emission as is the case for the AD model. The pa- 




rameters of the accretion flow are ?n out = 8 x 10"~ 4 , 
r tr = 50, S = 0.3 and s = 0.6. The parameters of 
the jet are rhj e t — 1.2 x 10“ 5 , p = 2.05, e e — 0.6 and 
e B = 10“ 4 (Pjet “ 4.8 x 10 42 erg s~ 3 , Pjet/Pboi ~ 23 and 
m ie t/m(3P 5 ) « 3 x 10” 5 ). 

We adopted two different transition radii in the fits 
shown in figure 5. The accretion flow r model with rtr = 
360 overpredicts the UV data (left panel of fig. 5). A 
smaller transition radius of v tr = 50 on the other hand 
does a better overall job of describing the optical-UV 
data. Such small transition radius seems also to be fa- 
vored by the modelling of the Fe Ka emission line (Young 
et al. 2007). We note however that with the available 
near-IR and optical data it is not possible to set a strong 
constraint on the properties of the truncated thin disk. 

M81 has been the subject of several modeling efforts 
(Quataert et al. 1999; Markoff et al. 2008; Yuan, Yu & 
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Ho 2009), Our goal here is not to discuss extensively this 
particular LINER since this was done elsewhere (Miller 
et al 2010). but rather put the general features of our 
AD and JD models for M81 in context with the results 
found by other authors. 

Quataert et ah (1999) did not consider the jet contribu- 
tion and favored an AD scenario for the origin of the SED 
(they used an “old” ADAF model assuming <5 — 0.01 and 
s = 0, i.e. no winds). Markoff et al. (2008); Yuan, Yu & 
Ho (2009) on the other hand were also able to explain the 
observations, but in the context of a JD model. Yuan, 
Yu & Ho (2009) also considered the simultaneous contri- 
bution of the ADAF and jet but restricted themselves to 
8 = 0.5 in their fits. Mainly for this reason. Yuan, Yu & 
Ho (2009) were unable to find an AD model for the SED, 
since as we mentioned before, ADAF X-ray spectra with 
such values of <5 are quite hard. 

It should be noted that Quataert et al. (1999) esti- 
mated a transition radius which is twice (r tr = 100) the 
value favored in the present work and a much higher 
accretion rate. The reason for this is that these authors 
previously fitted the SED of M81 using an outdated value 
of the black hole mass, which was fifteen times smaller 
and more uncertain than the value we use (estimated via 
spatially resolved gas and stellar kinematics, see EHF10). 
As a consequence, the multicolor blackbody radiated by 
their thin disk for the same accretion rate and transition 
radius that we use is considerably hotter and fainter. As 
a consequence, they needed to increase and m. 

4.5. NGC 3998 

The shortest wavelength (1750 A) UV data point in 
the SED of NGC 3998 presented in EHF10 is anoma- 
lously high because of variability (Devereux 2011). As 
discussed by Devereux this data point was obtained many 
years before all the other observations, when the source 
was much brighter. For this reason, we exclude this UV 
measurement from our analysis and plots. 

The left panel of Figure 6 shows an AD model in which 
the ADAF dominates the emission from the IR to the 
X-ray bands and explains well the available optieal-UV 
and X-ray data. As usual, the jet dominates the radio 
emission. The ADAF parameters are ?n ou t — 7.2 x 10" 3 , 
r out = 10 4 , (5 — 0.1 and s = 0.4. The jet parameters are 
m iet = 1.7 x 10" 6 , p = 2.2, e € - 0.01 and e B = KT 3 . 
For the resulting jet power from this model (i^ et = 9x 
10 42 erg s" 1 , rhj e t/rh(3Rs) ~ 6 x 10" 3 ), the jet kinetic 
energy is quite modest compared to the radiative output 
(Pjet/Tbol ~ 0.6). 

In the AD model above we adopted a large outer radius 
for the ADAF, but smaller radii are not ruled out by the 
data. For instance, we are able to obtain a reasonable 
fit to the SED with r tv = 500, m olit = 10" 3 and s = 
5 — 0.01, This model is consistent with the IR upper 
limits and accounts for the X-ray data. The transition 
radius cannot be much smaller than this value, otherwise 
the emission of the truncated thin disk would exceed the 
IR upper limits. Furthermore, as noted by Ptak et al. 
(2004), the lack of Fe Ka line emission also suggests that 
the value of r tr is not so small. 

The right panel of Fig. 6 shows a JD model in which 
the jet dominates the radio, optical-UV and the X-ray 
emission. Note that although the jet accounts quite well 
for the X-rays, it somewhat underestimates the optical- 


UV data. The parameters of the jet are rhj e t — 3.5 x 
10" 6 , p - 2.01, e e = 0.75 and e B = 3x 10" 5 (P jet - 
1.8 x 10 43 erg s" 1 , Pj e t/^boi ~ 1*3 and m ]et /rh(3Rs) « 
0.025). For illustration, we also show in the right panel 
an ADAF model computed with parameters such that 
its contribution in X-rays is weak (m out — 4 x 10" 4 , 
r tr = 100, S = 0.3 and s = 0.3). Note that, as for 
the cases discussed in the previous section, for values of 
8 > 0,3 the ADAF X-ray spectrum is harder than the 
data. 

The SED of this NGC 3998 was modelled before by 
Ptak et al. (2004). Here we confirm the main results of 
their work, namely that both AD and JD models are able 
to account w r ell for the broadband SED. 

5. FITS TO SPARSELY SAMPLED SEDS 

We present in this section the model fits to the 16 SEDs 
that were not discussed in the previous section. 

The SED fits are shown in Figures 7-12. In the sub- 
sections below, we briefly describe the details of the fits 
to the SED of each object. The model parameters are 
summarized in Table 2. 

5.1. NGC 0266 

The AD model accounts slightly better for the best-fit 
X-ray slope, but note the pronnounced uncertainty on 
the value of observed photon index. The JD model re- 
quires p < 2 which is below the usual range 2 < p < 3 
suggested by relativistic shock theory (Bednarz & Os- 
trowski 1998; Kirk et al. 2000). 

5.2. NGC 1097 

Note that even though the SED of NGC 1097 is as 
well sampled as many of the LLAGNs in Section 4, we 
chose to include it in this section rather than discussing 
it in more detail in §4 because its SED was already well- 
studied by Nemmen et al. (2006). 

We take into account here the HST UV spectrum which 
was not included in EHF10 because it contains contri- 
butions from an obscured nuclear starburst in addition 
to the AGN. This spectrum was obtained from Storchi- 
Bergmann et al. (2005) who studied and modeled it in 
detail. 

Our models do not take into account the contribution 
of the nuclear starburst. Therefore, the AD and JD mod- 
els do not reproduce the “UV bump” between « 1000 A 
and « 7000 A ( see Fig. 3 in Storchi-Bergmann et al. 
2005). 

The AD model displayed in Figure 7 corresponds to 
the model obtained by Nemmen et al. (2006). The JD 
model is also able to account for the whole SED, though 
it does not reproduce the slope of the X-ray spectrum 
as well as the AD model. Regarding the JD fit, given 
the available data there is no need to incorporate the 
contribution of the ADAF, Hence we do not include any 
ADAF model in the right panel of Fig. 7. 

5.3. NGC 1553 

The Bondi accretion rate was estimated bv Pellegrini 
(2005). 

Even though the JD model roughly accounts for the es- 
timated Lx, it fails to fit the slope of the X-ray spectrum 
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even with a small value of p . We take this as evidence 
that this source is unlikely to he JD. The estimated m ou t 
is consistent with the lower limit on rheondi obtained by 
Pellegrini (2005). 

5.4. NGC2681 

For this LLAGN, there are not enough optical obser- 
vations to fit the truncated thin disk model and estimate 
the transition radius. 

5.5. NGC 3169 

As was the case of NGC 2681, for this LLAGN there are 
not enough optical observations to fit the truncated thin 
disk model and estimate the transition radius. The JD 
model accounts well for the available data. The inferred 
extreme values of e e and e b imply that essentially all the 
energy in the post-shock region of the jet is carried by 
the particles. 

5.6. NGC 3226 

The SED of this LLAGN is well fitted by both an AD 
and JD type of models. As is the case of NGC 2681, there 
are no good optical band constraints on the emission of 
the truncated thin disk. 

5.7. NGC 3379 

This source has only two data points outside the X-ray 
band, one in the radio and the other in the optical, both 
of which correspond to upper limits. Therefore, there are 
few constraints for the accretion-jet model. The Bondi 
rate was estimated by David et al. (2005). 

Figure 9 shows that the AD model is able to reproduce 
the observed SED given the few constraints available, 
but there are no radio data to constrain the jet model 
in this case so we don't include it. The available data 
are not enough to constrain well the transition radius. 
The accretion rate required by the AD model is more 
than an order of magnitude higher than mBondi- One 
possible explanation for this result is that the accretion 
rate is enhanced by gas released by stars. This is in line 



with the findings of Soria et ah (2006a,b) for a sample of 
quiescent early- type galaxies. 

5.8, NGC 

The contribution of the truncated thin disk is required 
in order to account for the optical emission. The required 
transition radius in particular is r tr = 70. 

The JD model requires a value of p smaller than the 
usual values suggested by relativistic shock theory, as was 
the case of NGC 266. 

5.9. NGC 4261 

For this object, estimates of the Bondi rate, inclination 
angle and jet power are available (Gliozzi, Sambruna & 
Brandt 2003; Merloni & Heinz 2007). There are two in- 
dependent estimates of the jet pow r er. The first is based 
on energetic considerations and the VLBA observations 
of the pc-scale jet, which yield a lower limit to the jet 
power (Gliozzi, Sambruna & Brandt 2003). In the sec- 
ond method w r e use the 5 GHz specific luminosity and 
the Merloni & Heinz (2007) correlation obtaining a much 
higher jet power estimate. 

Note that the accretion rate that w r e obtain from the 
AD fit is somewhat smaller than raeondi- This AGN is 
likely to be strongly affected by extinction in the OUV 
band (EHF10). Hence, the OUV measurements probably 
do not capture the emission of the central engine. For 
this reason, it is not surprising that the models overpre- 
dict the OUV emission by almost one order of magnitude 
in Fig. 9. 

5.10. NGC 4278 

The radio emission of NGC 4278 was previously stud- 
ied in the context of ADAF models by Di Matteo, Carilli 
& Fabian (2001), who found that s > 0 is required in 
order for the radio emission predicted by ADAF models 
not to overestimate the radio flux. Di Matteo, Carilli 
& Fabian (2001) also estimated that the Bondi rate is 
mBondi ^ 6.001 - 0.01. Giroletti et al. (2005) found a 
two-sided radio structure for the jet using VLA data, 





log(i/L„ / erg s' 1 ) log(i >L V / erg s^ 1 ) \og(vL„ / erg s- 1 ) 


12 


Nemmen et al. 








Figure 7. SEDs and coupled accretion-jet models for NGC 0266, NGC 1097 and NGC 1553. The left panel shows the AD models for 
each object while the right panel displays the JD models. The dashed, dotted and dot-dashed lines correspond to the emission from the 
ADAF, truncated thin disk and jet, respectively. The solid line when present represents the sum of the emission from all components. 
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Figure 8. Same as Figure 7 for NGC 2681, NGC 3169 and NGC 3226. 
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Figure 9. Same as Figure 7 for NGC 3379, NGC 4143 and NGC 4261. 
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and estimated that the jet is oriented close to the line of 
sight (2° < i < 4°) and mildly relativistic (Tj ~ 1.5). In 
our jet models we therefore adopt i = 3° and Tj = 1.5. 

Both the AD and JD models require a truncated thin 
disk with a small transition radius (r tr ~ 30 ~~ 40) in 
order to account for the near-IR data. The values of the 
transition radius that we used in the models presented 
above (r tr = 30 — 40) are not unique and models with 
larger radii can also reproduce the SED. For instance, 
an accretion model with r tr = 100, rh Q ut = 4 x 10" 3 , 
5 == 0.1 and s = 0.77 is able to explain the optical and 
X-ray observations. 

The jet emission by itself under predicts the 1 mm radio 
observation, therefore the ADAF contribution is required 
even in the JD case. In the JD fit, a small value of s is 
favored by the observations as is the case of the JD fit 
for NGC 4143. 

5.11. NGC 4457 

Since there are only upper limits in the radio band, 
the jet power for this object was estimated using the 
observation at v = 1.5 x 10 10 Hz and the Merloni & 
Heinz (2007) correlation. 

5.12. NGC 4494 

Same as NGC 4457. 

5.13. NGC 4548 

Same as NGC 4457. Given the large uncertainty in the 
photon index, the X-ray spectrum does not provide good 
constraints for the models. 

5.14. NGC 4552 

The Bondi accretion rate was estimated by Merloni & 
Heinz (2007) (see also Allen et al. 2006) from the X-ray 
profiles of density and temperature. The kinetic power 
carried by the jet was estimated by Merloni & Heinz 
(2007) (see also Allen et al. 2006) from the energy de- 
posited in the X-ray cavities. 

Since the X-ray spectrum is quite soft, AD models are 
unable to account for the X-ray emission. The JD model 
in Fig. 11 is roughly consistent with the radio observa- 
tions and explains quite well the X-ray data, but overpre- 
dicts the optical data. The resulting jet power is roughly 
consistent with the value estimated by Allen et al. (2006); 
Merloni & Heinz (2007). 

5.15. NGC 4519 

NGC 4579 shares many characteristics with NGC 3031 
(M81) and NGC 1097: its nucleus features broad double- 
peaked Balmer emission lines (Barth et al. 2001) and a 
lack of the iron Ka line emission (Eracleous et al. 2002). 
Quataert et al. (1999) previously modeled the SED of 
NGC 4579 compiled by Ho (1999) using the “old” ADAF 
model (Narayan & Yi 1995). 

From the width of the broad Ha line, Barth et al. 
(2001) obtained a rough estimate of the inner radius of 
the line-emitting portion of the accretion disk r tr ~ 160. 
Hence, we adopt in our SED models r tr < 160 and 
i ~ 45°. 

As can be seen in Fig. 11, the optical continuum data 
do not require a small transition radius since it can be 
accounted by either the jet (JD fit) or ADAF (AD fit) 


emission. The truncated thin disk in our fits is some- 
what hotter than the thin disk modeled by Quataert et 
al. (1999), since their adopted black hole mass is out- 
dated and is ~ 10 smaller than the value we use in this 
work. The transition radius can be chosen to be some- 
what smaller if the accretion rate is lower and still be 
consistent with the data. 

We show together with the JD model an ADAF model 
which is the same as the one calculated for the AD fit 
except for the higher value of s. 

5.16. NGC 4756 

As was the case of NGC 1097, NGC 4143 and NGC 
4278, this LINER requires a small transition radius in 
order to explain the OUV data. 

6. THE DISTRIBUTION OF THE ACCRETION- JET 
PARAMETERS FOR LLAGNS IN LINERS 

Our detailed SED modeling allows us to place indepen- 
dent constraints on the mass accretion rate onto the black 
hole, the jet mass-loss rate and jet kinetic power for the 
LLAGNs in our sample. These quantities are of direct 
relevance to studies of the feeding and feedback of the 
underfed supermassive black holes in nearby galaxies. It 
is therefore prudent to study the resulting distributions 
of the central engine parameters from our models, 

6.1. Accretion rates 

In the previous section we showed that essentially two 
types of models (AD and JD) are able to account for the 
SEDs, depending on the combination of parameters that 
we choose. We found that in most of the JD fits, the con- 
tribution of an underlying ADAF is not required at all 
(cf. the JD models for' NGC 4594, M87 and NGC 4579). 
Since the estimate of the accretion rate relies on fitting 
the ADAF-thin disk model to the data, in those cases no 
robust rh constraint is available. We therefore show in 
the left and right panels of Figure 13 the histograms dis- 
playing the distribution of values of m ou t (dimensionless) 
and Afout (in units of Mq year" 1 ) respectively, consider- 
ing only the AD models. For one object (NGC 4552) we 
are unable to find a suitable AD model, therefore we do 
not include this object in the histogram. 

Most of the objects have accretion rates in the range 
rh « 0.003 — 0.01. The mean accretion rate for the sam- 
ple, in the case of the AD fits, is (logm) = —2.5 or 

/log^M/Me- vear" 1 ^ = —2 with standard deviations 

of 0.5 and 0.6, respectively. 

Not all gas supplied at the outer radius of the accre- 
tion flow ends up being accreted onto the black hole due 
to mass-loss via winds produced in the ADAF. From 
our derived parameters we can estimate the fraction of 
the gas that actually falls into the event horizon. Fig- 
ure 14 shows the histogram of values of M(r isco ) tak- 
ing into account the radial dependence of the density 
profile in the ADAF, wdiere r{ SCO = 3. We obtain 

^log ^Af(n sco )/M0 year" 1 ^ — ~3.3 with a standard 

deviation of 0.5. In other words, typically only a frac- 
tion of 14% of the available fuel supply reaches near the 
ISCO and presumably a similar fraction is accreted onto 
the black hole. 
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Figure 10. Same as Figure 7 for NGC 4278, NGC 4457 and NGC 4494. 
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Figure 11. Same as Figure 7 for NGC 4548, NGC 4552 and NGC 4579. 
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logM/M^a logM(M e yr- 1 ) 


Figure 13. Distribution of accretion rates at the outer radius of the ADAF obtained from our AD SED fits, a, left: Accretion rate in 
Eddington units, b, right: Accretion rate in physical units. 


6.2. Jet powers 

Especially relevant for studies of the feedback of mas- 
sive black holes in the local universe is the amount of 
energy that the jets deposit in their environment. Fig- 
ure 15 shows the distribution of jet kinetic powers ob- 
tained via both the AD and JD type of SED fits. As w r e 
noted above, we were not able to estimate a jet power for 
the AD model fit to NGC 4552. The average jet power 
/log (Pjet/erg s" 1 )) is 41.8 and 42.2 respectively for the 
AD and JD scenarios, each with a standard deviation of 
« 1 dex. Therefore, we can say that the typical jet power 
is ^ 10 42±1 erg s"" 1 . The distributions in Fig. 15 are in 
qualitative agreement with the distribution of jet powers 
estimated by Nagar, Falcke & Wilson (2005) (see their 
Fig. 6). 

Also relevant in studies of the kinetic AGN feedback 
is what fraction of the rest mass energy associated with 
the mass accretion rate is tapped by the central engine 
and converted into the jet kinetic power, which can be 


defined as r)$ et = jFj et /( Me 2 ). The choice of accretion 
rate in the definition of the “jet kinetic efficiency” is ar- 
bitrary and we could choose to use either the accretion 
rate M out at the outer radius of the ADAF or the one at 
the innermost stable circular orbit (ISCO), Mi sco (note 
that Mi SCO < Mout given the mass- loss through winds in 
the ADAF). The ISCO radius in this case corresponds to 
3 Rs, appropriate for a Schwarzschild black hole. Figure 
16 shows the distribution of the jet kinetic efficiency de- 
fined using Af 0 ut (top panel) and M\ sco (bottom panel) 
for the case of AD SED fits, since in the JD fits there 
are no reliable estimates of the accretion rates. The av- 
erage values of Iogr/j et are -2.9 (M = M ou t) and -2.1 
(M = Mi SCO ) (both these averages have an uncertainty 
of « 0.7 dex). Clearly, the latter efficiency is higher in 
order to be consistent with the smaller ISCO accretion 
rates. 

We are also able to quantify the fraction of the mass 
accretion rate that escapes the attraction of the black 
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logM^} (m® y f 1 ) 


Figure 14. Distribution of accretion rates at the radius corre- 
sponding to the ISCO of a Schwarzschild black hole, estimated 
from the AD fits. 



Figure 16. Distribution of jet kinetic efficiencies Pj et /(Mc 2 ) con- 
sidering M = Mi sco and M = M 0 n t for the AD models. 
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Figure 15. Distribution of jet powers. 

hole and is channeled into the jet. Figure 17 shows 
the distribution of values of Mj et / Misco for the AD 
fits. As it turns out, the average value of the ratio is 

^log = —2.9, i.e. on average 0.4% of the 

material that reaches the ISCO is channeled in the jet. 
The corresponding uncertainty in this estimate is 0.6 dex. 



6.3. Bolometric luminosities and radiative efficiencies 

One possible way of estimating the bolometric lumi- 
nosities for the LLAGNs using the observed SEDs is to 
integrate the SEDs directly, not taking into account up- 
per limits (especially in the IR) and assuming a suitable 
interpolation. For instance, this procedure was carried 
out by EHF1G who also computed a bolometric correc- 
tion for the X-ray luminosity. 

In our case, we have physically motivated fits for 
the 21 SEDs. Therefore, for each LINER in our sam- 
ple, we calculated Lbd by integrating the total SED. 
Carrying this out using the SEDs obtained according 



n 

£ 


log M&IM^ 

Figure 17. Distribution of Mj e t/Misc<v 
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to the AD scenario we obtain the average value of 
(log (Lboi/ or g s" 1 )) ~ 41.5 with a standard deviation 
of « 1 dex. For the SED fits computed according to 
the JD scenario we find a corresponding average of 41.3 
with a similar standard deviation, which gives the same 
average as the one computed using the bolometric lu- 
minosities estimated by EHF10. Figure 18 displays the 
distribution of values of Lboi resulting from the models 
compared with the one estimated by EIIFIO, 

Defining the X-ray bolometric correction as kx ^ 
Lboi/Lx, we calculated the average bolometric correc- 
tion (log/cx) resulting from our models as 1.45 and 1.35 
according to the AD and JD scenarios respectively, with 
standard deviations of 0.36 dex (AD) and 0.46 dex (JD). 
The corresponding geometric means of kx are 28 (AD) 
and 22 (JD), the mean values are 43 (AD) and 47 (JD) 
and the medians are 22 (AD) and 14 (JD). Our models 
result in values of kx which are somewhat lower than 
those estimated by EHFIO. 
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Figure 18. Distribution of boiometric luminosities inferred from 
the models compared with the ones estimated by EHF10. 

Using our SED models we can calculate what is the 
typical fraction of the accreted matter that is converted 
into the combined radiation of the accretion flow and 
the jet, i.e. the radiative efficiency r] Vdb( $ = L^ 0 i/(M 0 Ut c 2 ) 
(note that the accretion rate in the definition of this ef- 
ficiency is the one at the outer radius of the ADAF), 
We obtained the average value (log7? ra d) = —3.25 with a 
standard deviation of 0.75 dex. Hence, the typical radia- 
tive efficiency is 5 x 10" 4 . 

We estimated how the jet kinetic power compares to 
the total radiative output of the central engine, i.e. we 
calculate the ratio P^ t /Lbo\ where P }et is obtained from 
the jet model and Lboi from integrating the modeled 
SEDs. We obtain (log (Pj et /Tboi)) ^ 0-2 for the AD-type 
of fits and 0.9 for the JD models (there is a « 0.5 dex 
uncertainty in both values). The reason for the smaller 
jet powers in the AD fits compared to the JD ones is 
that the mass-loss rates required in the former fits are 
systematically smaller than in the ones in which the jet 
has a higher degree of contribution in the X-rays. Based 
on our results and given the degeneracy in the fits, we 
can say that the typical kinetic to radiative power ratio 
of LLAGNs is « 5 (i.e. the average of the value from 
the AD and JD fits). This is in agreement with the find- 
ings of Nagar, Falcke & Wilson (2005) that the accretion 
power output in LLAGNs is dominated by the jet power. 
Figure 19 shows the distribution of values of Pj et /Lboi 
resulting from the models. 

7. THE AVERAGE SED 

It is of interest to compute the average SEDs resulting 
separately from the AD and JD scenarios. These SEDs 
are useful for different purposes. First and foremost, 
since the observed SEDs are the main observables that 
we use to derive the central engine parameters via our 
fitting, as a consistency check the average model SEDs 
should reproduce the average SED obtained directly from 
the observed ones. 

Secondly, SEDs obtained from averaging observed data 
points are obviously limited by the observed bands. Out- 
side ’the observed bands different authors adopt ad-hoc 



Figure 19. Distribution of the ratio of jet power to boiometric 
luminosity. 

interpolations between the data, usually a linear interpo- 
lation in log-log space (Ho 1999; EHF10), which might 
not reflect the actual physical processes involved in the 
emission. The average from the AD and JD scenarios 
that we obtained can hence serve as guidelines - with a 
physical justification - to the typical shape of the SEDs 
in the bands which have not been constrained yet. 

Motivated by the reasons above, we show in Fig. 20a 
the average SEDs computed separately for the AD and 
JD scenarios. We first normalized the individual SEDs 
to the same X-ray luminosity of 10 40 erg s'" 1 in the band 
2-10 keV. This value is approximately the average X-ray 
luminosity for the LINERs in our sample. After normal- 
izing the SEDs, for each frequency bin we computed the 
geometric mean of the model SEDs corresponding to the 
particular model class as (log^L^). Following EHF10 we 
choose to compute the average of the logarithm of the lu- 
minosities instead of using the values of vL u themselves 
in order to reduce the effect of outliers in the resulting 
average. 

The different bumps in the average model SEDs reflect 
the different physical processes that operate in the flows. 
In the average AD model SED, the different bands are 
dominated by different radiative processes in the flow as 
described below: 

1 m — 1 mm: jet optically thick synchrotron emission. 

1 mm — « 100/im: the first bump and usually the 
strongest one in the average AD SED is due to the 
ADAF synchrotron emission. 

« 10/mi — « 1/mi: the second bump is due to the ther- 
mal emission from the truncated thin disk. 

« Igm “0.1 keV: there is a weak UV bump due to 
the first inverse Compton scattering of synchrotron 
photons in the ADAF. 

« 0.1 keV — 100 keV: the last bump in X-rays is dom- 
inated by the second inverse Compton scattering of 
synchrotron photons in the ADAF. 
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The average JD SED is more simple and has only one 
hump in the IE between a few x 10 pm - a few xl/rni 
which is due to the truncated thin disk. The rest of 
the emission is due to synchrotron photons from the jet. 
Overall, the spectral shapes of the average JD and AD 
SEDs are quite similar in the radio and X-rays above 
1 keV. In between these bands the shapes are slightly 
dissimilar. The bumps in the average SEDs are much 
less pronounced than the bumps in the individual SEDs. 
This occurs because the bumps in the individual SEDs 
do not peak at the same place and when the SEDs are 
averaged these bumps are smoothed out. For this reason, 
the average model SEDs will not resemble any one of the 
individual SEDs. 

Figure 20a also displays the average data points com- 
puted from the observed SEDs in a similar way by EHF10 
where the error bars represent the uncertainty in the 
emission due to the uncertainty affecting the amount 
of extinction correction involved. The shaded region 
around the average best-fit X-ray power-law represents 
the standard deviation in the value of the LINER pho- 
ton index. In order to illustrate the wide diversity of the 
individual SEDs, we show in Figure 20b the la scatter 
affecting the model SEDs (where we show only the AD 
average SED for simplicity) and the observed ones. 

As expected, the average model SEDs agree well with 
the observed constraints. There are some details that 
are worth mentioning. The shape of the X-ray spectrum 
of the average JD SED is slightly softer than the corre- 
sponding shape of the average AD SED. Both are within 
the la uncertainty in the photon index of the average ob- 
served SED. In the OUV, even though the model SEDs 
agree with the observed constraints they are quite dif- 
ferent from each other. For instance, the red bump is 
stronger in the average JD SED. The average JD SED 
predicts a lower level of UV flux. In the radio band, the 
model SEDs are quite similar to each other. 

Figure 21 shows the average AD and JD SEDs com- 
pared to the average ones of radio-loud and radio-quiet 
quasars computed by Shang et al. (2011). The average 
quasar SEDs computed by Shang et al. (2011) are very 
similar to the ones of Elvis (1994) but the former in- 
clude more detailed features, and are based on more re- 
cent data obtained with improved instrumentation. As 
in Fig. 20a, the SEDs were normalized such that they 
all have the same X-ray luminosity in the 2-10 keV band 
of 10 40 erg s" 1 . The UV excess in the quasar SEDs (the 
big blue bump) is clearly apparent in comparison to the 
LLAGN ones. It is also interesting that for v > 10 1 ' Hz 
the average AD, average JD and the radio- loud quasar 
SEDs are quite similar. 

8. DISCUSSION: THE NATUEE OF X-RAY 
EMISSION 

In our AD models, the X-rays are produced predom- 
inantly by inverse Compton scattering by the ADAF of 
seed synchrotron photons produced in the accretion flow 
itself; in contrast, in our JD models, the X-rays are dom- 
inated by the optically thin synchrotron emission in the 
base of the jet. In this regard, the Chandra X-ray spec- 
trum is one of the most important constraints for the 
origin of high-energy photons in LLAGNs. 

There are only two objects for which the AD or the 
JD can be ruled out, within the observational uncertain- 


ties and the explored model parameter space: NGC 1553 
(§5.3) and NGC 4552 (§5.14). For NGC 1553, the JD 
model does not reproduce the hardness of the X-ray spec- 
trum whereas for NGC 4552 the AD model is not able 
to account for the corresponding X-ray softness. For the 
remaining objects, by inspecting the SED fits we can see 
that both the AD and the JD models are able to repro- 
duce with varying degrees of success the whole SED from 
radio to X-rays, in particular the X-ray portion. 

In order to quantify which model systematically fits 
better the broadband SEDs we use as a goodness-of-fit 
parameter a modified y 2 . We define this modified y 2 
simply as 

X — (l/obs “ y mo del) (2) 

i 

where y = \ogvL v and the sum is carried over the ob- 
served data points. We sample the X-ray continuum us- 
ing 3 data points from 0.5 keV to 10 keV and we do not 
take into account upper limit data points. The values of 
y 2 for each model are displayed in Table 2. 

Considering the 20 objects for which we were able to 
obtain both JD and AD fits, we find that y 2 (JD) < 
X 2 {AD) for 9 objects and y 2 {JD) > y 2 (AD) for the 
remaining 8. In other words, for 9 objects (10 including 
NGC 4552 for which we did not find an appropriate AD 
fit) the JD results in better SED fits as opposed to 8 for 
which the AD model is a better fit. The results from 
this goodness-of-fit test show that 52% of the SEDs in 
our sample are better fitted by the JD model whereas 
38% are better reproduced by the AD scenario. For the 
3 objects with y 2 (JD)/y 2 (AD) « 1 (NGC 1097, NGC 
3998 and NGC 4374) it is hard to tell - based on the 
broadband SED alone - which model provides a better 
fit. 

We should point out a number of caveats involved when 
considering which scenario better accounts for the nature 
of the X-ray emission. From an observational point of 
view, even though we are carrying the largest systematic 
analysis of the SEDs of LLAGNs in LINEEs done up to 
date, the size of our sample is limited. We need to extend 
our modeling to a larger number of LLAGNs in order to 
increase the statistics and draw further conclusions. Sec- 
ondly, for many sources we only have a few data points 
available to fit, so clearly for these objects we need to 
obtain more measurements to better constrain the mod- 
els. In addition, in our modeling we find that one of 
the main observables that help to pinpoint the dominant 
source of X-rays is the hardness of the X-ray spectrum. 
ADAF models tend to predict harder X-ray spectra while 
jet models predict softer ones. In most of the LLAGNs 
in our sample there is a significant observational uncer- 
tainty in the X-ray hardness derived from the Chandra 
observations which make it difficult to determine which 
whether the AD or JD is favored based on the X-ray 
spectrum. 

From a theoretical perspective, there are considerable 
theoretical uncertainties involved in the the ADAF-jet 
models. For instance, the values of the model param- 
eters that affect the microphysics and dynamics of the 
ADAF and jet are not sufficiently constrained. In par- 
ticular, S and p have a large impact on the hardness of 
the X-ray spectrum predicted by ADAF /jet models but 
there is a considerable uncertainty in the range of pos- 
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Figure 20. a, left: The average SEDs (geometric mean) computed separately for the AD and JD models (dashed and dotted lines 
respectively). The data points correspond to the geometric mean computed by EHF10. b, right: la scatter around the average (model 
and observed) SEDs illustrating the diversity of individual SEDs. The solid line shows the average AD SED and the shaded region 
corresponds to the standard deviation from the AD models. The points correspond to the mean computed by EHF10 and the error bars 
show the scatter in the measurements. In the OUV, the filled circles correspond to measurements without any reddening correction whereas 
the open circles correspond to the maximal extinction correction. 



Figure 21. The average JD and AD model SEDs compared to 
the average radio-loud and radio-quiet quasar SEDs computed by 
Shang et al. (2011). 


sible values of these parameters despite recent progress 
(e.g., Sharma et al, 2007; Medvedev 2006). Considering 
the overall goodness-of-fit of the AD and JD models, it 
is important to point out that our jet model is basically 
a phenomenological one. Uncertainties affecting for ex- 
ample the jet model could translate in larger values of 
A 2 . For instance, when we find x 2 (J D) < x 2 (AD) 1 that 
could be possibly because the JD model is less well con- 
strained than the AD one, i.e, we have more freedom in 
the jet fitting compared to the ADAF model. 

Taking into account these caveats, it is still premature 
to favor one model over the other based only on the SED. 


The nature of the X-ray emission in LLAGNs has been 
debated in the last few years by several authors favoring 
in some cases the AD or JD models based on the analysis 
of individual sources (e.g., Falcke & Markoff 2000; Yuan, 
Markoff & Falcke 2002; Yuan et al. 2002; Yuan, Quataert 
& Narayan 2003; Wu, Yuan & Cao 2007; Markoff et al. 
2008; Miller et al. 2010) or in a statistical sense based 
on the fundamental plane of black hole activity (Mer- 
loni, Heinz & Di Matteo 2003; Falcke, Kording & Markoff 
2004; Yuan & Cui 2005; Yuan, Yu & Ho 2009; Plotkin 
et al. 2011). We suggest different routes to clarify this 
issue in the future. 

The AD and JD models should predict different charac- 
teristic radio and X-ray variability timescales (e.g., Ptak 
et al. 1998), hence we should be able to test which com- 
ponent is dominant in X-rays by carrying out a simul- 
taneous monitoring of the variability of radio and X- 
rays. By comparing the variability pattern predicted by 
jet/ ADAF models we could pinpoint the X-ray dominant 
component. This promising strategy is quite similar to 
what the observational campaigns carried out for M81 
(Markoff et al. 2008; Miller et al. 2010) and the radio- 
quiet Seyfert 1 NGC 4051 (King et al 2011) and could be 
applied to many other LLAGNs. We note that in the case 
of M81, Miller et al. (2010) favor a scenario in which the 
X-rays originate in the transition region between the thin 
disk and the ADAF; NGC 4051 is a different case since 
it is a relatively bright AGN accreting at the ~ 0.05LEdd 
level (Peterson et al. 2004) and its jet-disk coupling is 
quite distinct compared to LLAGNs. 

The Fermi Gamma-ray Space Telescope can be quite 
helpful in this regard. LLAGNs are potential sources of 
7 -rays (Takami 2011) and in fact one of the source in 
this sample (NGC 4486) has already been detected by 
the Fermi LAT (Abdo et al. 2009). By using the mode! 
parameters that we derived in our radio- to- X-rays fits we 
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should be able to predict the 7- ray spectrum (Mahade- 
van, Narayan & Krolik 1997; Takami 2011) and compare 
with future Fermi detections. In this way we should be 
able to compare AD/JD predictions and through 7-ray 
observations place further constraints on the production 
site responsible for the high-energy emission and the jet- 
disk connection in LLAGNs. 

Observations in the mm and sub-rum are also very 
helpful because they constrain the ADAF synchrotron 
emission and hence also the X-ray emission. This fol- 
lows because the synchrotron photons in the ADAF are 
inverse-Compton-seattered to X-rays. 

Finally, further refinements and a better understand- 
ing of the fundamental plane of black hole activity (Mer- 
loni, Heinz & Di Matteo 2003; Falcke, Kording & Markoff 
2004; Giiltekin et al. 2009; Yuan, Yu & Ho 2009; Plotkin 
et al. 2011) should make the fundamental plane a better 
tool for constraining the radiative processes shaping the 
radio and X-ray emission in sub-Eddington black hole 
source and by consequence LLAGNs. 

9. SUMMARY 

We performed detailed modeling of the broadband - 
radio-to-X-rays - spectral energy distributions of a sam- 
ple of 21 low-luminosity AGNs in LINERs selected from 
EHF10. With this exploratory modeling, our goal is to 
constrain the general properties of the central engines of 
the dominant population of AGNs at z ~ 0. Our coupled 
accretion-jet model consists of an accretion flow which is 
radiatively inefficient in the inner parts and becomes a 
thin disk outside a certain transition radius. The rela- 
tivistic jet is modeled in the framework of the internal 
shock scenario. 

We demonstrated that there are two classes of mod- 
els that can explain the majority of the observed SEDs. 
We call the first one AD which stands for “ADAF- 
dominated” since the ADAF dominates most of the 
broadband emission, particularly the X-rays. In the sec- 
ond class of models, the jet component dominates the 
majority of the continuum emission and for this reason 
we call this scenario JD as in “jet-dominated” . The SEDs 
predicted by the AD models have a more complex shape 
than the JD ones, given the richer variety of radiative 
processes involved in ADAFs. For both scenarios though, 
the radio band is almost alw r ays dominated by the syn- 
chrotron emission from the jet and the near IR to optical 
band is dominated by the truncated thin disk thermal 
emission. 

From our exploratory modeling of the SEDs we are able 
to constrain important parameters that characterize the 
central engines of LLAGNs, as summarized below. 

Mass accretion rates: based on the AD models, the 
values of the accretion rates supplied to the accre- 
tion flow lie in the range 4 x 10~ 4 — 0.02 in Ed- 
dington units with a mean (logm) = —2.5 or alter- 
natively /log ^M/M.0 year” 1 ^ } = — 2 (standard 

deviations of « 0.5 dex respectively). Of that gas 
supply, typically only ~ 10% is accreted onto the 
black hole with the remaining gas being lost due to 
outflows. 

Jet powers: the typical jet power resulting from the AD 


and JD scenarios is ~ 10 42 erg s 1 ranging from 
~ 10 40 to ~ 10 44 erg s" 1 . 

Jet production efficiencies: The typical efficiency 
with which the rest mass energy associated with 
gas supplied to the accretion flow is converted into 
jet kinetic power is estimated to be ~ 0.001 (range 
~ 10“ 4 — 0.01). The efficiency is ten times larger if 
we consider the gas that reaches a few gravitational 
radii only. 

Fraction of material ejected in the jet: On aver- 
age, a fraction of 0.4% of the gas that reaches a few 
gravitational radii is chanelled in the relativistic 
jet (range 0.01 — 1%). 

Radiative efficiencies: The typical efficiency of con- 
version of rest mass energy associated with gas sup- 
plied to the accretion flow into disk-fjet radiation is 
« 5x 10 -4 (dispersion of 0.8 dex), producing an av- 
erage bolometric luminosity of a few x!0 41 erg s" 1 . 

Kinetic to radiative output ratio: LLAGNs typi- 
cally release « 5 times more kinetic power in the 
jets than the total energy radiated away. 

Transition radii between thin disk and ADAF: A 

truncated thin disk is required to fit the available 
mid and near-IR data for only 4 sources - NGC 
1097 (see also Nemmen et al. 2006), NGC 4143, 
NGC 4278 and NGC 4736 - with the resulting 
transition radii in the range 30 — 22 5R$. 

The values that we derived from our SED models above 
provide useful indicators of the feeding and feedback 
properties of LLAGNs in LINERs, and by extension of 
the whole LLAGN population. 

Even though we used the physical scenario that is fa- 
vored to account for the combined set of observational 
properties of LLAGNs (Ho 2008; Narayan & McClintock 
2008), the inferred values of the parameters are model de- 
pendent. There are theoretical and observational sources 
of uncertainty in the values of the derived parameters. 
For instance, there are degeneracies between the effect of 
the model parameters on the resulting SEDs (Quataert & 
Narayan 1999b; Nemmen et al. 2006). The uncertainty in 
the OUV - due to uncertain degree of extinction correc- 
tion - and the X-ray data also contribute to uncertainty 
in the estimated parameters. Furthermore, our sample is 
currently limited to 21 objects. In spite of these caveats, 
the estimates we obtained for the properties of central 
engines of LINER AGNs should be illustrative of the gen- 
eral behavior of the class. 

We can draw a link between the supermassive black 
holes in LINERs and the quiescent black hole in our 
galaxy, Sgr A*. By modelling the SED of Sgr A* with 
ADAF models, an accretion rate of m ou t ~ 10 5 was es- 
timated (Yuan, Quataert & Narayan 2003; Yuan, Shen & 
Zhi-Qiang 2006). This value is two orders of magnitude 
below the typical accretion rate of LINERs that we esti- 
mated. Therefore, we could say that in order to Sgr A* 
“light up” and become a LINER-like LLAGN, it needs 
to accrete at a rate ~ 100 times higher than the present 
one. 
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The individual and average SED models can be useful 
for a different number of applications. They provide a li- 
brary of templates for the compact emission of LL AGNs, 
which can be useful for studies of the emission of other 
astrophysical systems in the nucleus, he., stellar popula- 
tions, dust and PAH features (e.g., Storchi-Bergmann et 
ah 2005; Mason et ah 2007). Furthermore, the models 
provide the expected emission on bands that were not 
vet observed in many LLAGNs (e.g.. mid and far IR. far 
UV). 
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0.002 
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Notes: 

(a) Black hole mass estimated from the fundamental plane of black holes (Merloni, Heinz k Di Matteo 2003), using the X-ray and radio luminosities. 

(b) Observed jet power estimated using the radio data and the Merloni k Heinz (2007) correlation. 

(c) For NGC 1097, V = 10 (for consistency with the work of Nemmen et al. 2006). 

(d) 5 GHz luminosity estimated using the Merloni, Heinz k Di Matteo (2003) correlation. 

(e) The lower limit on PP* >S was derived by Gliozzi, Sambruna k Brandt (2003) and the upper limit results from using the radio data and the Merloni k Heinz (2007) correlation. 
The Bondi rate and inclination angle is from Gliozzi, Sambruna k Brandt (2003). 

References: 1. Pellegrini (2005), 2. David et al. (2005), 3. Gliozzi, Sambruna k Brandt (2003), 4. Di Matteo, Carilli k Fabian (2001), 5. Giroletti et al. (2005), 6. Allen et al. 
(2006), 7. Merloni k Heinz (2007), 8. Biretta, Sparks k Macchetto (1999), 9. Di Matteo et al. (2003), 10. Barth et al. (2001) 
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Table 2 

Model parameters resulting from the SED fits. The meaning of the parameters is described in Section 3. 
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Table 3 

Model parameters resulting from the SED fits (continued). 
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